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ABSTRACT 

We study the polycyclic aromatic hydrocarbons (PAH) bands, ionic emission lines, and mid- 
infrared continuum properties, in a sample of 171 emission line galaxies taken from the literature 
plus 15 new active galactic nucleus (AGN) Spitzer spectra. We normalize the spectra at A = 23/im 
and grouped them according to the type of nuclear activity. The continuum shape steeply rises 
for longer wavelengths and can be fitted with a warm blackbody distribution of T ~ 150-300K. 
The brightest PAH spectral bands (6.2, 7.7, 8.6, 11.3, and 12.7 /im) and the forbidden emission 
lines of [Sin] 34.8 /im, [Aril] 6.9 /im, [Sill] 18.7 and 33.4 /im were detected in all the starbursts 
and in ^80% of the Seyfert 2. Taking under consideration only the PAH bands at 7.7/zm, 11.3/xm, 
and 12.7/im we find that they are present in ^80% of the Seyfert 1, while only half of this type of 
activity show the 6.2/im and 8.6 PAH bands. The observed intensity ratios for neutral and 
ionized PAHs (6.2/im/7.7/im x 11.3/j,m/7.7/j,m) were compared to theoretical intensity ratios, 
showing that AGNs have higher ionization fraction and larger PAH molecules (> 180 carbon 
atoms) than SB galaxies. The ratio between the ionized (7.7/zm) and the neutral PAH bands 
(8.6/xm and 11.3/mi) are distributed over different ranges for AGNs and SB galaxies, suggesting 
that these ratios could depend on the ionization fraction, as well as on the hardness of the 
radiation field. The ratio between the 7.7/zm and 11.3^m bands is nearly constant with the 
increase of [Ne in] 15.5/im/[Ne n] 12.8/zm, indicating that the fraction of ionized to neutral PAH 
bands does not depend on the hardness of the radiation field. The equivalent width of both PAH 
features show the same dependence (strongly decreasing) with [Nein]/[Nen], suggesting that the 
PAH molecules, emitting either ionized (7.7/im) or neutral (11.3/xm) bands, may be destroyed 
with the increase of the hardness of the radiation field. 



Subject headings: galaxies: Seyfert 
spectroscopic 



laxies: starburst - infrared: ISM - ISM: molecules - techniques: 



1. Introduction 

The Mid-infrared (Mid-IR) spectra of galaxies 
either with active nucleus (AGN) and/or Star- 
burst (SB) show emission features attributed 
to Polycyclic Aromatic Hydrocarbons (PAH), 
which can be considered to be originated in very 
small amorphous carbon dust grains or very large 
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carbon-rich ring molecules (e.g. iPuget fc Leger 



19891: iDraine fe Leel Il984t braine fe Li I l200ll ) 



The most prominent, well-known PAH emissions, 
are the 6.2/i m 7.7jnm, 8.6jnin, lPg ^ m and 12.7Mm 



bands 



Roche et al. 



Weedman et al. 



, ,19911 
2005; 



Genzel et al. 



Buchanan et al. 



1998; 

2006). Besides, the Mid-IR spectral region 
also presents prominent forbidden ionic emis- 
sion lines, such as [Neil] 12.8/im, [Nev] 14.3/im, 
[Ne in] 15.5^m, [S in] 18.7ami and 33 48/nm, [O iv] 25.89ami 
and [ Si n] 34.8/xm (e.g. lSturm et al. 1120021: IBernard-Salas et al. 
120091) . 

The nature of the ionizing source can be as- 
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sessed by computing the intensity ratios be- 
tween forbidden lines such as [N eil] 12.8/xm , 
[Nem] 15.5^m and [Nev] 14.3/zm (e.g. lVoit Ill992t ). 
In addition, the detection of [Ne v] at 14.3/im and 
24.2^m or [Ne vi] at 7.6/im, due to ions of high 
ionization potential, is an evidence of a hard radi- 



ation field associated wi th an AGN (jGenzel et al. 
19981 : ISturm et al.ll200d ). 



Several diagnostic diagrams based on the line 
intensity ratios between the brightest ionic lines 
and the PAH features have been proposed to 
classify the Mid-IR spectra of galaxies accord- 
ing to the degree of activity (AGN, low-ionization 
nuclear emission-line r egion - LINER, and SB 



galaxies). For instance, Genzel et al. ( 19981) 



us- 
ing data from the Short Wavelength Spectrom- 
eter (SWS) on board of the Infrared Space Ob- 
servatory (ISO), propose the diagnostic diagram, 
[O iv] 25.9/zm/[Nen] 12.8/xm versus the strength 
of the 7.7/Ltm PAH band to separate star-forming 
galaxies from AGNs. 

From the study of the Mid-IR spectra of 
quiescent spirals, Galactic Hn regions, plan- 
etary n ebulae, and p hotodissociation regions 
fPDRs'l . [Galliano I (120061) found that the ratios be- 
tween the 6.2/im/8.6/im and 7.7/im/8.6/im PAH 
bands do not vary significantly, while the ratios 
6.2^m/11.3/xm and 7. 7/im/11.3/^m varies over one 
order of magnitude. They concluded that the 
properties of the PAHs for different types of galax- 
ies are very similar and that variations of the PAH 
band ratios are due to the fraction of ionized to 
neutral PAH molecules. Thus, this ratio could 
be used to analyze the physical conditions of the 
molecular gas in the regio ns where the emission 
is originated. In addition, iGordon et al. (2008) 
found that the equivalent width (EW) of PAH fea- 
tures, observed in Hn regions of the M101 galaxy 
are correlated with the ionization index, and does 
not depend on the Hn regions metallicity. 

According to ISmith et al. (2007) the ionized 
to neutral fraction of PAHs can be estimated 
using PAH 7.7/xm/11.3/xm ratio, which is sensi- 
tive to the type of nuclear activity, decreasing 
from SB to Seyfert. These authors, suggest that 
this behavior is due to the selective destruction 
of the 7.7/im PAH molecule with the increasing 
hardness of the radiation field (similar results are 
found bvlO'Dowd et al. 1 120091: Hunt et al. Il2010t 
Kaneda et al. 20081) . Such results are supported 



by laboratory experiments, which show that cer- 
tain PAHs are effectively destroyed by individual 
UV and X-ray photons, and cannot survive within 
a few kiloparsecs from an active nucleus, unless 
they are shielded from the AGN's X-ray emission 
by absorbing material ( Voit Ill992h . 

It is worth mentioning that the above results, 
concerning the ionization fraction of PAH bands 
in AGNs are carried out with small samples of 
this type of object. Therefore, we use the largest 
sample of AGNs and SB galaxies available up to 
date in Spitzer public archive, to study in detail 
the effect of the hardness of the radiation field on 
the PAH bands emission, in these kind of galax- 
ies. This paper is structured as follows: in Sect. 
[2] we describe the sample selection, data reduction 
process and emission line measurements. In Sect. 
[3] we present results and discussions on the shape 
of the continuum spectra, emission line frequency 
and diagnostic diagrams of PAH band ratios and 
emission lines ratios. Sect. [4] we discuss the be- 
havior of the ionization fraction and EW of PAH 
bands with the hardness of radiation field. Sum- 
mary and conclucion are presented in Sect. [5] 

2. The Sample 

We have analyzed a well defined sample of 
186 galaxies which have Spitzer spectra available 
from public archive. B eing 83 AG Ns taken from 
Gallimore et~aTl (|201dh . 22 SB from llBrandl et al ' 
20061). 59 Hn and LINER from (jsmith et al. 
20071 ) and 7 Hn re gions of M101 galaxy from 



(|Gordon et al~ll2008h . plus 15 AGNs Spitzer spec- 
tra published here for the first time. Details on the 
reduction process of the 15 new objects are given 
below. 

2.1. Data Reduction 

The 15 new spectra of AGN galaxies (4 Sy 1 
and 11 Sy 2) were obtained from the Spitzer public 
archive. Observations wer e performed with the In- 
frared Spectrograph (IRS; iHouck et al. us- 
ing the Short-Low (SL) and Long-Low (LL) mod- 
ules, covering the interval between 5fim and 36/im 
with a resolving power of 64 - 128. The raw data 
wereprocessed using the Spitzer pipeline version 
17.2_|. The nuclear spectra were extracted using 



1 Available at: http: / /ssc. spitzer. caltech.edu/spitzerdataarchives/ 
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the Spitzer IRS Custom Extractor - SPICEo De- 
fault point-source extractions were selected for all 
objects. Only a few sources show a noticeable step 
in flux between the overlapping LL and LS mod- 
ule spectra, which results from the extended cir- 
cumnuclear emission contribution to the size SL 
module slit of 3 .7 x 57 and LL module slit of 
10 .7 x 168 . Although the observed regions of 
almost all of the galaxies are large in both mod- 
ules the AGN light dominates the spectra. Fig. 
[1] shows the emission spatial profile of NGC3786. 
Note that the luminosity falls off to almost zero for 
« 3 pixels from the center of the nucleus. For this 
galaxy, 3 pixels correspond to « 1 kpc, adopting 
a radial velocity of 2678 km sec -1 and a Hubble 
constant of 74 km sec _1 Mpc _1 . 

The general parameters for the objects listed 
in Table CD were taken from NASA/IPAC Extra- 
galactic Database (NED). All the spectra were 
corrected for rcdshift. The rest-frame spectra were 
grouped according to their activity, and are shown 
in Fig. 1 

2.2. Emission Line Measurements 

For the 15 ne w AGNs liste d in th e Tableland 

~ (|2006l ). low resolu- 



SB taken from Brandletal 



tion IRS spectra have been decomp osed using the 
PAHFi iH code, described in detail bv lSmith et al 



( 20071 ). In short, they assume that the IRS spectra 
are composed by dust continuum, starlight, promi- 
nent emission lines, individual and blended PAH 
emission bands, and that the light is attenuated 
by extinction due to silicate grains. 

The default input parameters of pahfit are: (i) 
the infrared emission from a blackbody with T= 
5000 K; (ii) different weights for thermal dust con- 
tinuum components represented by the tempera- 
tures 35, 40, 50, 90, 135, 200, 300K; (iii) spectral 
line features, including the pure rotational lines of 
the molecular hydrogen; (iv) dust features are rep- 
resented by individual and blended Drude profiles, 
which are theoretical frequency profiles for a clas- 
sical damped harmonic oscillator and are the best 
choice to model PAH emission. For more details 



sec 



Smith et al. I (|2007l ) 



With the above constraints we were able to iso- 



2 Available at: http:/ /ssc. spitzer. caltech.edu/postbcd/spice. html. 
3 Source code and documentation for PAHFIT are available in 
http://tir.astro.utoledo.edu/jdsmith/research/pahfit.php 



late the PAH features and measure their fluxes, 
as well as those of the ionized gas emission lines. 
The fluxes and the EW are listed in Tables H to [3] 
We checked the consistency of the results with the 
liner code ( Pogge fc Owen 19931 ) that integrates 
the flux of a Gaussian function fitted to the line, 
obtaining similar results in both procedures. The 
observed IRS spectra of Mrk334, a Sy 1 galaxy, is 
shown in Fig. 01 where we also plot the detailed 
components of the pahfit spectral decomposition 
model. Clearly, the model reproduces well the ob- 
served spectrum. 

3. Results and Discussion 
3.1. New spectra 

The main goal of this section is to charac- 
terize the continuum shapes and spectral fea- 
tures observed in the 15 new objects listed in 
Table [TJ with respect to the type of nuclear ac- 
tivity (Sy 1 and Sy 2). Therefore, we normal- 
ize to unity the continuum of all spectra at A 
= 23 /im and grouped them according to nu- 
clear activity. For each AGN class, data were 
sorted according to spectral shape, from bluest 
to reddest, top to bottom of Fig. [3J respectively. 
One can clearly see that the majority of spec- 
tra have a continuum shape increasing for longer 
wavelengths, in agreem e nt with results found b y 



Deo et al. 



Weedman et al 

~iE009l) 



(2005); Buchanan et al. I (2006) 



( 20091) 



Wu et al. 



. . , _ .. .Burtscher et al. . , _ 
(|2009t) ; lBaum et all (|20Rf 7 This continuum range 
can be fitted with a warm blackbody distribu- 
tion of T ~ 150 — 300 K. Interestingly, a simi- 
lar trend in the continuum of AGNs is also ob- 



served in the near-infr ared region (jRiffel et al 
20061 : iRiffel et al. II2009I) . 



We compare the spectroscopi c properties of the 



15 new objects with those of the lBuchanan et al. 



( 2006h groups. We found that the spectra of the 
Sy 2 galaxies: NGC3786, NGC5728, NGC7682, 
Mrk471, Mrk609, Mrk622, Mrk883, Mrkl066, 
Mrk883 and the Sy 1 Mrk334 and NGC4748, have 
the same characteristic as their group 1 (strong 
PAH emission and red continuum, see Fig [2J . On 
the other hand the Sy 2, NGC 1275, NGC2622, 
Mrk3 and Sy 1 Mrk 478 can be classified as group 
3 (power-law continuum and weak PAH features), 
while only the Sy 1 Mrk279 galaxy show a broken 
continuum and silicate emission at 9.8/im (group 
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2). I n summary, for our and Buchanan et al. I 
( 20061 ) samples the Mid-IR spectra for most of Sy 
2 galaxies have same spectroscopic characteristics 
as those of group 1 (> 50%). 

3.2. PAH and Emission Line Frequency 

Fig. [3] shows a histogram of the commonest 
emission features present in each type of activity 
(note that we have used the 186 galaxies spectra). 
It is clear from this figure that the PAH features 
are present in almost all the Sy 1, Sy 2 and SB 
galaxies, as well as the H2 molecules at 9.6/xm and 
17/Um. In addition, the most frequent ionic lines 
are: [Sin] 34.8/xm, [Aril] 6.9/im, [S in] 18.7 and 
33.4/xm, [Siv] 10.5/xm, [Ncm] 15.53/xm and [Nev] 
14.3/im. 

All studied SB galaxies show PAH bands and 
the following forbidden emission lines: [Si 11] 34.8/im, 
[Aril] 6.9^m, [S in] 18.7 and 33.4/im. In addition, 
nearly 80% show molecular H2 lines. Similarly, 
>80% of the Sy 2 show PAH features, and the 
same ionic lines detected in SB galaxies. However, 
we have found that PAH lines in Sy 1 galaxies be- 
have differently from the other types of nuclear 
activity. The PAH bands at 7.7/im, 11.3/im and 
12.7^m are present in 80% of the galaxies, but 
only 50% of them show the PAH lines at 6.2/im 
and 8.6/im. Nevertheless, the high ionization lines, 
such as [Nev] 14.3^m and [Oiv] 25.8/im are more 
common in Sy 1. 

3.3. Diagnostic Diagram of PAH bands: 
grain sizes and ionization fraction 

Numerical studies of the vibrational energy dis- 
tribution of PAH molecules show that the relative 
strengths of the 3.3/zm, 6.2/xm, 7.7/mi, 8.6/im. and 
11.3^m PAH features depen d on grain size an d 



on the charging conditions ([Draine fc Li I [2001) 



Small PAHs radiate strongly at 6.2/im and 7.7/im, 
while la rge PAHs emit mos t ly at longer wave- 
lengths (ISchutte et al. Ill993t iDraine k Li 1 12007 : 
Tielend 20081) . Ano ther important conclusion of 
Draine fc Li I (120011) is that neutral PAHs have 
higher values of 11.3^m/7.7/xm ratio, while ionized 
PAHs have lower ratios. On the other hand, the 
6.2/xm/7.7/xm ratio decreases with the increasing 
number of carb on atoms of the PA H molecules, 
Fig. 16 from IDraine fc Li I (|2001l ) shows how 
the relative strengths of these emission vary de- 



pending on the size and c harge state of PAH 
molecules. lO'Dowd et al. (2003) used this dia- 
gram to analyze these PAH ratios for a sample 
formed by galaxies dominated by star formation, 
galaxies with significant stellar and AGN compo- 
nent and AGN-dominated galaxies. They found 
a weak trend in the direction of constant ioniza- 
tion fraction with changing grain size. In addi- 
tion, they suggest that the presence of an AGN 
component is correlated with the reduction in the 
7.7/xm/11.3/um ratio indicating that smaller PAH 
grains (number of carbon atoms <180) are de- 
stroyed in AGNs. 

We present in Fig. [5] the theoretical in- 
tensity ratios, for neutral and ionized PAHs 



(6.2/xm /7.7/im x 11.3^m/7.7/im). from lDraine fc Li 
( 200ll ). The authors calculated these quanti- 
ties assuming grai ns illuminated either by the 



Mathis et al.l (|1983l ) spectrum or a blackbody with 



T= 3 x 10 4 K varying the number of carbon atoms 
in each molecule. We overplot in this diagram 
the observed PAH emission line ratios of our sam- 
ple. Clearly, the observed ratios are located be- 
tween neutral and ionized PAH theoretical lines, 



in agreement with previou s studies ([Draine fc Li 



200ll : IO'Dowd et al. Il2009h 



The above results suggest that the observed 
PAH bands in both, AGNs and low ionization 
objects (SB and LINERS), are formed by an ap- 
propriate mixture of PAH molecules with differ- 
ent sizes and adequate neutral to ionized frac- 
tion. High ionization objects (AGNs) present 
a tendency for molecules > 180 carbon atoms 
(6.2/xm/7.7/xm < 0.2), on the other hand most 
of low ionization objects are biased towards 
molecules with less than 180 carbon atoms. The 
latter are typically located in the region 0.2 
<6.2/xm/7.7/xm < 0.4 (Fig. [5]). Furthermore, 
Seyfert galaxies appear to be located near to re- 
gion populated by ionized PAH molecules, while 
the majority of SB, Hn, and LINER are near to 
the line representing the neutral PAHs. Moreover, 
these results are consistent with a picture that 
AGN have higher ionization fraction and larger 
PAH molecules than SB galaxies. 

Another way to assess the ionization fraction 
of the PAH molecules is using the following diag- 
nostic diagram 6.2/im/11.3/im x 7.7jnm/ 1 1.3jnm . 



Such diagram was used by Galliano ( 2006) , 



Galliano et al~l (|2008l ) and lO'Dowd et al 



2009), 
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they demonstrate that there is a good correlation 
between both ratios. However, we would like to 
emphasize the fact that both samples were mostly 
composed by SB galaxies. In addition, these ratios 
are co ntrolled by the fraction of ionized to neutral 



PAHs (jAllamandola et al.lll999h 



In Fig [5] (top) we show these diagram for our 
sample. Clearly, if we only consider the SB galax- 
ies our results are consistent with the correla- 
tion described above. By including the AGNs 
(Sy 1 and Sy 2), an unexpected large dispersion 
is observed for both type of activity. Alterna- 
tively, we also test the diagram 6.2/im/8.6/im 
x 7.7/xm/8.6/xm involving ionized (7.7/im and 
6.2/xm) and neutral (8.6/im) PAH molecules. This 
diagram is shown in Figj6] (bottom), where a sim- 
ilar dispersion is observed. Interestingly, in both 
diagnostic diagrams the ratios in the vertical axes 
tend to separate the AGNs (Sy 1 and Sy 2) from 
SB galaxies. This is even clearer in the diagram 
involving the 7.7/xm/8.6/im ratio, leading us to 
the conclusion that 7.7/.tm/8.6/.tm > 6 for AGNs 
(dotted line) and smaller for low ionization ob- 
jects. 

Synthetic spectra in the 6.2/im - 9.0/mi region, 
for large symmetric PAHs catio ns, anions and neu- 
tral, h ave been computed by iBauschlicher et aL 
(|2008l ). The trends, in the band position and 
intensite, is a function of the molecule size, 
charge and geometry. Large PAH cations and 
anions (> 110 carbon atmos) produce prominente 
bands at 7.7/im and 8.5/im, the intensite ratio 
7.7/zm/8.6/zm increase with the PAH sizes. Thus, 
we interpret the separation in activity types (SB, 
Sy 1 and Sy 2), along the vertical axis in both dia- 
grams of Fig. [51 as due to the fact that in average 
the emitting PAH molecules in Seyferts are larger 
than those of SB. In addition, we suggest that 
the high values for the ratios 7.7/im/8.6/im and 
7.7/im/11.3^m is due to the fact that almost all 
the molecules are ionized. In summary, we point 
out that both diagrams (Figs. [5] and [6]) clearly 
shows that the ionization fraction and molecules 
size increase from SB to AGNs. However, it is 
worth mentioning that, no separation between 
Seyfert type is observed. 

Although, the four ratios depends on the PAHs 
ionization fractions, ratios involving 6.2^m band 
have the same range of variation for both types 
of nuclear activity. This can be associated with 



the fact that th e 6.2/xm PAH emission is due to 
small molecules ( Bauschlicher et al. 2008l h which 
probably are destroyed by the hard radiation field 
of Seyfert galaxies. 



Diagnostic Diagram for Emission Line 
Objects 



3.4. 



As shown in Sect. 13.21 (Fig. [3]), forbidden ionic 
emission lines are observed in SB and active galax- 
ies. In the first case, the gas is ionized by massive 
hot stars, while in the latter the ionization is due 
to a non-thermal continuum. Moreover, the emis- 
sion cloud in A GNs are located in the narrow line 
regions (NLR) (|Osterbrock fc Ferland Il2006l ). 

The Mid-IR emissi on from dusty NLR o f AGNs 



has been explored bv iGroves et al.l (j2006l ). Their 



models give a direct connection between the dust 
emission and the line emission from the photoion- 
ized gas. These authors also show that the emis- 
sion line ratios [Nev]/[Neii] x [Nem]/[Neii] can 
be u sed to separate AGNs from SBs. In addi- 
tion, iThornley et al. I ( 2000l ) using star burst evo- 
lutionary models, which incorporate new stellar 
atmosphere models for massive stars, found that 
[Ne in] / [Ne 11] is sensitive to hardness of the radi- 
ation field. 

In Fig. [7] we show the diagnostic diagram 
[Nev]14.3/im/[Neii]12.8/xm x [Nelil]15.5/xm/ 
[Ne n]12.8/im for our sample. We also overplot 
the NLR models for differen t ionization and pres - 
sure parameters taken from IGroves et al.l ( 2006 ) . 
Clearly, all the observations have values lower than 
those predicted by the models. In this diagram, 
AGNs are closely correlated, while the SB galaxies 
show a large scattering. 

The [Ne v] 14.3/zm has a high ionization poten- 
tial of 97.1 eV and is often detected in spectra of 



AGN s (e.g. ISturm et al~l 120021 IWeedman et~al 



20051) however is rarely detected in SB galaxie s 



Brandl et ahl 120061 iBernard-Salas et al. M2009I) 



The lack of correlation for SB galaxies in this dia- 
gram, thus, may be due to a strong contamination 
of [Ncv] with PAH 14.19/im emission or to the 
absence of the former line in the SB galaxies spec- 
tra. In addition, we call the attention to the fact 
that the deblending of these lines is very difficult 
in the low resolution Spitzer spectra. 

It is also clear in this figure that Sy 1 galax- 
ies have [Nem] / [Neil] > 0, while SB and Sy 2 
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have this ratio < (solid line in Fig. [7]). This 
suggests that both line ratios are sensitive to the 
hardnes s of the radiation field, confirming the re- 



sult of iThornlev et al. I (|2000t) . The fact that 



Sy 2 galaxies are spread over the lower left of the 
Fig. [7] can be associated to iGroves et al. (l2006h 



predictions, which states that starburst contribu- 
tion would move the points in that direction. 

It is worth mentioning that there are 3 SB 
galaxies located in the top-right side of Fig. [7J 
(NGC1097, NGC4676 an d NGC520). The former, 
was re classified as Sy 1 bv lStorchi-Bergmann et aL~l 
(jl997h . while the two latter are strongly interact- 
ing source s and may have the thes e emission lines 
enhanced ( Read fc Ponman]|l998l ). Thus we sug- 
gest that further, more accurate studies (i.e. with 
higher spectral resolution spectra), are necessary 
to better explain the position of the SB galaxies 
in the diagram. 

4. Behavior of the Ionization Fraction and 
EW of PAH Molecules with the Hard- 
ness of Radiation Field 



Smith et al. ( 20071 ) have studied the correla- 



tion between the intensity ratio of the strongest 
PAH bands 7.7/im/11.3/im with the hardness of 
the radiation field indicator, [Neni] / [Neil] . They 
have found that galaxies with Hn regions or star- 
burst optical spectra exhibit a nearly constant 
7.7/im/11.3/xm across the whole range of radia- 
tion hardness, including low-metallicity galaxies, 
which have [Nem] /[Neil] < 2. In contrast, the 
AGNs are located below the Hn-type galaxies 
and their PAH ratio falls rapidly with increas- 
ing i^cUatic^i_^eld.^_This result was confirmed 
bv lO'Dowd et al~ ( 20091 ) through the analysis of 
PAH emission bands of a sample of 92 typical star 
forming galaxies. They interpreted the decrease in 
the PAHs band ratio for AGNs as due to a selective 
destruction of the ionized 7.7/^m PAH molecules. 

The variation of 7.7/xm/11.3/xm PAH bands, 
observed in our sample (186 objects), with respect 
to the hardness of the radiation field, is shown in 
Fig. [H It is clear from this figure, that the fraction 
of ionized to neutral PAH bands does not depend 
on the hardness of the radiation field in contrast 



wit h the previous r e sults of ISmith et al. I ([2007) 



PAH bands would have the same dependence on 
the radiation field. 

Besides, this diagram show a tendency in the 
sense that SB galaxies have [Nem] /[Neil] < 0.4, 
while for Sy 1 and Sy 2 [Nem] /[Ncn] > 0.4. In 
addition, Sy 1 are biased towards location on the 
upper right corner. Some of the Sy 2 galaxies are 
located in the SB region, which suggest that these 
galaxies may have a circumnuclear starburst. No 
clear separation between Sy 1 and Sy 2 is observed 
for [Nem] /[Ncn] > 0.4. 

In order to test the dependence of the PAH 
bands strength with the radiation field, we com- 
pare in Fig. [3] the EW of the 7.7/zm and 11.3/im 
PAH bands with [Nem] /[Neil] hardness indica- 
tor. It can be seen that EW of both PAH bands 
are constant with the radiation field indicator for 
[Ne III] /[Ne n] < 0.8, and fall rapidly with the in- 
creasing radiation field. Galaxies, powered by star 
formation, have a nearly constant EW across the 
full range of [Nem] /[Nell] . While for Sy 1 the 
EWs of the PAH bands decrease rapidly, for the 
Sy 2 the EW values are distributed either in the SB 
or Sy 1 region. Their location in the diagram may 
reflect the stellar or AGNs dominant component 
contribution to the spectra (i.e. in some case the 
Sy 2 may be dominated by a non-thermal compo- 
nent, and in other by circumnuclear massive star 
formation). In order to have a more quantitative 
analysis of Fig. |H1 we have computed an expo- 
nential regression taking into account all objects 
(AGNs+SBs, dotted line) and another for a sub- 
sample formed only by AGNs (solid line). For the 
latter we have found a correlation coefficient of ~ 
-0.7 while for the former we have not found any 
correlation (~ -0.2). 

We conclude this section by arguing that both 
PAH molecules, either ionized (7.7/im) or neutral 
(11.3/zm), may be destroyed with increasing hard- 
ness of the radiation field. We emphasize that sim - 



ilar results were found by iGenzel et all (|l998h . 



which compares the PAH strength with the high 
to low excitation line ratios for a sample of ULIRS 
(ultra lu minous infrar e d gal axies) and AGNs. In 
addition, Baum et al. ( 20101) has shown that the 
diagram EW PAH 6.2/xm x [Ne v] 14.3/xm/[Ne n] 
12.8^m separate AGNs from SB galaxies. 



and lO'Dowd et al. I ([20091 ) . We interpret this be- 
havior as due to the fact that the strengths of both 
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5. Summary and Conclusions 

We have analyzed the emission line properties 
of a sample of 186 Spitzer galaxy spectra. 171 
objects were taken from the literature and we 
have added 15 unpublished AGNs Spitzer spec- 
tra, taken from the public archive. The continuum 
properties of the 15 new spectra were also studied. 
Our main results can be summarized as follows. 

• The normalized Mid-IR continuum spectra 
(at A = 23/im), of the 15 new AGNs, steeply 
rise for redder wavelengths and can be fitted 
with a warm blackbody distribution of T ~ 
150-300K. 

• The brightest PAH spectral bands (6.2, 
7.7, 8.6, 11.3, and 12.7/um) and the forbid- 
den emission lines of [Sin] 34.8 /xm, [Aril] 
6.9 /im, [S in] 18.7 and 33.4 /im were de- 
tected in all the SB and in ~80% of the 
Sy 2. Taking under consideration only the 
PAH bands at 7.7/im, 11.3/im, and 12.7/im 
we find they are present in ~80% of the 
Seyfert 1, while only half of this type of 
activity show the 6.2 / um and 8.6 /zm PAH 
bands. 

• Comparison between observed line inten- 
sity ratios of neutral and ionized PAHs 
(6.2fim/7.7fim x 11.3/im/7.7/im) with the- 
oretical models suggests that these ratios, in 
Sy 1, Sy2 and SB galaxies, are produced by 
an appropriate mixture of PAH molecules, 
with different sizes and adequate neutral to 
ionized fraction. 

• The PAH emission lines observed in Sy 1 and 
Sy 2 may be due to PAH molecules with > 
180 carbon atoms, while for most of the low 
ionization objects (SB and LINERs) it is bi- 
ased towards molecules with < 180 carbon 
atoms. In addition, Seyfert galaxies appear 
to be located near to region populated by 
ionized PAH molecules, while the majority 
of SB, Hn, and LINER populate the region 
of neutral PAHs. These results are consis- 
tent with a picture where Sy 1 and Sy 2 have 
higher ionization fraction and larger PAH 
molecules than SB galaxies. 

• We have investigated the PAH line inten- 
sity ratio diagrams 6.2^m/11.3^m x 7.7/im/ 



11.3/xm and 6.2^m/8.6^m x 7.7/im/8.6/im. 
The separation in activity types (SB, Sy 1 
and Sy2), along the vertical axis in both dia- 
grams, is interpreted as due to the fact that 
in average the emitting PAH molecules in 
Seyferts are larger than those of SB, with 
ionization fraction and emitting molecules 
sizes increasing from SB to AGNs. No sepa- 
ration between Seyfert type is observed. 

• Diagnostic diagram [Ne v]14.3/im/[Ncii]12.8^im 
x [Nein]15.5/^m/[Neli]12.8/xm shows that 
AGNs are correlated, while SB galaxies show 

a large scattering. Moreover, for most of the 
Sy 1 [Ncin] /[Neil] > 0, while for SB and 
Sy 2s [Ne in] /[Neil] < 0, indicating that 
both line ratios can be used as an indicator 
of radiation field hardness. 

• The ratio between the 7.7/Ltm and 11.3//m 
PAH bands is nearly constant with the in- 
crease of [Ne III] 15. 5^im/[Nen] 12.8^m, indi- 
cating that the fraction of ionized to neutral 
PAH bands does not depend on the hardness 
of the radiation field. 

• The equivalent width of both PAH features 
(7.7^im and 11.3/xm) show the same depen- 
dence with the [Ne in] /[Nell] ratio. In the 
case of Sy 1, they are nearly constant for 
[Nem] /[Neil] < 0.8, and fall sharply with 
increasing radiation field. For SB galaxies 
the PAH EW is constant with the increase of 
[Ne in] /[Neil]. Sy 2s are distributed either 
in the SB or Sy 1 regions. These results sug- 
gest that the PAH molecules emitting either 
ionized (7.7^m) or neutral (11.3//m) bands, 
may be destroyed with increasing hardness 
of the radiation field. 
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Fig. 2. — Calibrated spectra of the Sy 1 (top) and Sy 2 (bottom) galaxies. All spectra were ordered according 
to their shapes from a steeper spectrum (top) to a natter one (bottom). Some emission lines arc identified. 
All spectra were normalized at 23fim. 
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Fig. 3. — Frequency histogram of most common Mid-IR emission lines (in percent). SB, Sy 2 and Sy 1 are 
represented by empty and shaded bars, and solid lines, respectively. 



Iteration 102, Chi-Sq= 1 .1 3e + 03, Tau_9_7 = 0.030E Iter ation 137, Chi-Sq= 5.72e+03, Tau_9_7= 0.Q2 41 




Wavelength (furi) Wavelength (^m) 



Fig. 4. — Detailed decomposition of Mrk334 from 5 to 37/im using pahfit code. Solid black lines represent 
thermal dust continuum components. Measurements are represented by squares, uncertainties are plotted 
as vertical error-bars, which are smaller than the symbol size. Dotted black line indicates mixed extinction 
components. Continuous lines (gray) represent the best fit model with individual PAH features and emission 
lines, respectively. 
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Fig. 5. — Diagnostic diagram of 6.2/im/7.7/im x 1 1.3Atm/7.7Atm PAH ratios. The long dash line represents 
theoretical intensities of neutral to ionized PAHs ( Draine fc Li I [2001 ) . Dashed-dotted line shows the po- 
sition of molecules forme d by 180 carbon atm os. Empty triangles are Hii and L INER objects taken from 



Smith et al. I (|2007l ) and iGordon et al. 



_ (120081), full triangles ar e SB t aken from iBrandl et all (|2006l ). full 
and empty circles are Sy 1 and Sy 2 of our and iGallimore etUl fl2010h sample. Error bars are shown only 
for the 15 new AGNs, similar error were found for the remaing objects. 
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Fig. 6. — Diagnostic diagram of the relative strengths of the neutral (8.6/mi and 11.3/im) compared to ionized 
(6.2^m and 7.7/xm) PAH bands. In the top, 6. 2 /xm /11.3^ m x 7.7Mm/11.3Atm PA Hs. Empty triangles are Hn 



and LINER objects take n from Smith et al. I ( 20071 ) and Gordon et al. (120081). full triangle s are S B taken 



from lBrandl et al.1 ( 20061 ). full and empty circles are Sy 1 and Sy 2 of our and lGallimore et all ( 20101 ) sample. 
Error bars are shown only for the 15 new AGNs, similar error were found for the remaing objects. Dotted 
line separate high and low ionization objects. 
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Fig. 7. — Diagnostic diagram of [Nev] 14.3^m/[Neii] 12.8^m x [Nem] 15.5^m/[Neii] 12.8^m ratios 
direction of increasing ionization parameter for dusty NLR models with the total pressure of: P tot /k 
(dot c urve), ~ 10 7 (dashed curve) and ~ 10 s (solid curve) . Full triangles are SB taken from iBrandl et al 



The 
10 6 



( 2006f ). full and empty circles are Sy 1 and Sy 2 of our and iGallimore et all (|2010h sample. Error bars are 
shown only for the 15 new AGNs, similar error were found for the remaing objects. 
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Fig. 8. — Diagnostic diagram between an indicator of the hardness of the radiation field, 
[Nem] 15.5/[Ne li] 12.8mhi x 7.7/xm/ 11.3 pn PAH emission fea tures. Empty triangles are Hn a nd LINER ob- 
j ects t aken from lSmith et al." ] (|2007l i and bordon et al. I d2008h. full triangles arc SB taken from Brandl et al 
(|2006l ) , full and empty circles are Sy 1 and Sy 2 of our and iGallimore et al. I (|2010h sample. Error bars are 
shown only for the 15 new AGNs, similar error were found for the remaing objects. 
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[Nelll] 15.5pm/[NeII] 12.8|am 



Fig. 9. — Diagnostic diagram between an indicator of the hardness of the radiation field, 
[Nem] 15.5/ [Nell] 12.8£im ve rsus EW 7.7/xm (top ) and EW 11.3//m (bottom) PAH emission features. Full 



triangles are SB take n from iBrandl et al.l (|2006l ). full and empty circles are Sy 1 and Sy 2 of our and 
iGallimore et"aD (l2010h sample. Error bars are shown only for the 15 new AGNs, similar error were found for 
the remaing objects. The lines represent an exponential regression fitted for the total sample (dotted line) 
and for AGNs sub-sample (solid line). 
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Table 2: Aromatic Emission Line Strengths (10 W m ) and Equivalent Widths (/xm) 



Mrk270 2.88±0.33 0.035 7.19±0.S1 0-095 15. 6± 1.82 0.220 
Mrk334 28.5±0.30 0.548 9.92±0.67 0.183 107±2.88 1.881 



0.09±0.26 0.001 0.84± 0.12 0.013 2.94± 0.36 0.049 - 0.15±0.10 0.002 6.62±0.89 0.154 

20.2±0.24 0.351 1.12± 0.09 0.019 28. 3± 0.35 0.479 8.05± 0.19 0.131 16. 7± 0.39 0.263 18. 5± 0.83 0.229 



Sy 2 



Mrk3 0.77±0.31 0.006 4.56±0.63 0.037 38.2 ± 2.98 0.315 

Mrk471 2.67±0.12 0.820 3.98±0.32 1.030 6.78 ±1.15 1.415 

Mrk609 16.9±0.22 1.020 8.57 ±0.49 0.485 58.0 ±1.74 3.115 4.62± 0.20 0.245 

Mrk622 5.01 ±0.19 0.516 2.26±0.42 0.194 14.6 ± 1.65 0.986 

Mrk883 4.62±0.14 1.100 4.68±0.35 0.920 10.7 ± 1.18 1.630 

Mrkl066 47.0±0.28 0.757 15.5±0.58 0.263 199 ±2.76 3.567 18.6±0.37 0.340 

NGC1275 10.8±0.25 0.097 18±0.48 0.139 - 

NGC2622 2.75±0.14 0.233 4.33±0.37 0.339 - 

NGC5728 10.6±0.32 0.379 4.11±0.64 0.167 65.1 ± 3.00 3.050 

NGC7682 40.4±0.16 1.840 11.5±0.35 0.456 160 ±1.48 5.200 



I. 14±0.35 
2.15±0.09 

II. 7±0.16 
1.71±0.14 
1.33±0.10 
31.9±0.26 



8.09±0.31 
29.7±0.12 



1.22± 0.17 
0.05± 0.04 
0.53± 0.07 
0.37±0.07 
0.01±0.04 
1.29± 0.11 
10. 6± 0.10 
0.32± 0.05 
0.04± 0.21 
1.07± 0.05 



3.28± 0.55 
3.30±0.22 
17. 7± 0.32 
6.45± 0.28 
3.76± 0.22 
52. 1± 0.43 
17. 4± 0.33 
2.97±0.21 
21. 4± 0.65 
40. 9± 0.19 



1.61± 0.37 
0.56± 0.19 
5.76± 0.23 
2.40± 0.18 
0.83± 0.18 
15. 3± 0.23 
2.90± 0.22 
0.78± 0.14 
3.97± 0.35 
10. 1± 0.12 



1.45± 0.16 
1.67±0.29 
10. 4± 0.48 
4.41± 0.41 
1.37± 0.34 
30.9±0.45 

2.03± 0.13 
9.89± 0.78 
24.0±0.26 



39.2±1.57 
5.18±0.26 
16. 5± 0.66 
2.65±0.49 
2.53± 0.30 
43. 0± 1.41 
53.9±0.98 
2.42± 0.46 
32. 6± 1.14 
39.8±0.55 



4.88±0.37 0.245 
11.3±0.16 0.328 



0.009 
0.382 
0.617 
0.093 
0.161 
0.592 



0.424 
0.828 



0.009 
0.007 
0.027 
0.015 
0.001 
0.020 
0.043 
0.015 
0.001 
0.024 



0.022 
0.474 
0.903 
0.253 
0.324 
0.713 
0.070 
0.135 
0.783 
0.913 



0.009 
0.078 
0.282 
0.089 
0.068 
0.181 
0.011 
0.034 
0.117 
0.221 



0.008 
0.226 
0.474 
0.154 
0.107 
0.328 

0.086 
0.258 
0.522 



0.143 
0.359 
0.469 
0.044 
0.111 
0.325 
0.199 
0.0472 
0.471 
0.641 



Table 3: Atomic Emission Line Strengths (10 16 W m 2 ) and Equivalent Widths (/im). 

Name [Am 6.0;/m Ew [ Ar ; i ; 8 '. j >/<m Ew |\ 1 . 5 /i m Ew [Nc 1 1 1 2 . 8 ji m Ew [NenjJjXSjTm Ew > in lfs. 7/< m Ew O |\ [25 . 9jj in Ew [IVi; 2 r ) .9pm Ew i S : : : 3:t. L/i.m Ew 



Mrk279 0.70±0.16 

Mrk334 1.36±0.09 

Mrk478 0.47 ±0.03 

NGC4748 0.45 ±0.21 



0.014 0.05± 0.10 0.001 0.60±0.05 0.014 0.82± 0.05 0.022 1.26 ± 0.07 0.039 0.64 ±0.10 

0.036 0.06±0.05 0.001 0.63±0.03 0.016 5.32± 0.03 0.122 2.58±0.04 0.050 2.91±0.06 

0.012 0.06 ± 0.06 0.002 0.38±0.04 0.016 0.59 ± 0.04 0.032 0.52 ±0.01 0.032 0.40±0.01 

0.030 0.28 ± 0.12 0.018 2.34±0.09 0.144 1.14 ± 0.06 0.065 3.78 ±0.12 0.123 1.83±0.15 



0.024 1.10±0.05 0.067 

0.050 

0.027 0.75 ±0.02 0.066 

0.058 6.81 ±0.13 0.269 



0.87±0.13 0.088 

7.07±0.14 0.129 

0.01 ±0.06 0.001 

3.42 ±0.29 0.167 



0.49±0.06 0.008 



3.82±0.14 
- ±- 
2.69 ± 0.24 



Sy 2 



Mrk3 

Mrk471 

Mrk609 

Mrk622 

Mrk883 

Mrkl066 

NGC1275 

NGC3786 

NGC2622 

NGC5728 

NGC7682 



3.31±0.19 
0.43±0.06 
1.02±0.08 
0.40±0.08 
0.55±0.06 
3.75±0.12 
1.89±0.11 
0.56±0.06 
0.62±0.05 
0.95±0.10 
1.34±0.05 



0.040 2.63± 0.12 0.034 6.26±0.08 0.072 



8.75± 0.08 
0.32± 0.02 
2.22± 0.03 
0.70± 0.03 
1.31±0.02 
9.60± 0.03 
3.37± 0.05 
1.49± 0.02 
0.50± 0.02 
2.26± 0.06 
5.10± 0.02 



20.9 ± 0.12 
0.19±0.01 
1.11±0.04 
0.73±0.03 
0.96±0.03 
4.67 ± 0.09 
4.25 ± 0.14 
1.10±0.02 
0.91±0.03 
7.69±0.06 
3.79±0.04 



6.70 ± 0.11 
0.13±0.02 
1.19±0.05 
0.39±0.07 
0.71±0.03 
4.80±0.12 

1.77 ± 0.11 
0.06±0.03 
0.40±0.04 
3.85±0.08 
4.60±0.05 



18.2±0.10 
0.36±0.02 
0.82±0.03 
0.73±0.04 
0.91±0.02 
3.51±0.22 

1.31±0.19 
0.86±0.06 
13.7±0.09 
4.42±0.04 



-1.25±0.20 0.054 11.7±0.20 0.150 



0.153 
0.083 
0.046 
0.147 
0.096 
0.019 
0.031 
0.069 
0.061 
0.073 



0.11±0.02 
0.83±0.07 



0.11±0.04 
0.42±0.09 



0.019 
0.023 



0.009 
0.035 



0.12±0.01 
0.29±0.02 
0.34±0.02 
0.26±0.02 
1.61±0.04 
2.78±0.05 
0.51±0.02 
0.44±0.02 
2.84±0.10 
0.86±0.02 



0.028 
0.022 
0.021 
0.036 
0.038 
0.017 
0.031 
0.031 
0.197 
0.029 



0.073 
0.064 
0.144 
0.036 
0.150 
0.146 
0.019 
0.087 
0.031 
0.084 
0.164 



0.111 

0.020 
0.047 
0.019 
0.069 
0.055 
0.023 
0.042 
0.025 
0.162 
0.109 



0.037 
0.013 
0.048 
0.009 
0.043 
0.050 
0.009 
0.024 
0.012 
0.081 
0.093 



0.159 
0.049 
0.034 
0.020 
0.058 
0.031 

0.066 
0.036 
0.244 
0.077 



0.25± 0.29 0.002 



0.04 ±0.01 
0.23±0.05 



0.020 
0.009 



1.74±0.03 
1.59±0.08 
0.80±0.09 
0.82±0.05 
3.72±0.22 
0.07±0.18 
1.33±0.06 
0.59±0.05 
6.29±0.12 
6.03±0.09 



0.285 
0.070 
0.028 
0.052 
0.035 
0.001 
0.076 
0.036 
0.094 
0.102 



2.03±0.08 
3.63±0.13 
3.51±0.14 
2.58±0.09 
11.2±0.23 
10.2±0.23 
1.70±0.05 
2.56±0.08 
9.87±0.17 
12.2±0.11 



0.331 
0.159 
0.128 
0.161 
0.108 
0.096 
0.096 
0.165 
0.145 
0.203 



